Pairing state with a time-reversal symmetry breaking in FeAs-based superconductors 



Wei-Cheng Lee, 1 Shou-Cheng Zhang, 2 and Congjun Wu 1 

1 Department of Physics, University of California, San Diego, CA 92093 
2 Department of Physics, McCullough Building, Stanford University, Stanford CA 94305-4045 

We investigate the competition between the extended s±-wave and d x 2_ y 2-wa,ve pairing order 
parameters in the iron-based superconductors. Because of the frustrating pairing interactions among 
the electron and the hole fermi pockets, a time-reversal symmetry breaking s + id pairing state 
could be favored. We analyze this pairing state within the Ginzburg-Landau theory, and explore 
the experimental consequences. In such a state, spatial inhomogeneity induces supercurrent near a 
non-magnetic impurity and the corners of a square sample. The resonance mode between the s± 
and d a .2_ !/ 2-wave order parameters can be detected through the Bi 9 -Raman spectroscopy. 
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The discovery of the iron-based new superconductors 
with high critical temperatures has attracted a great deal 
of attention [H, [3, [3, [J H, H 0]. The symmetry struc- 
ture of Cooper pairs is one of the central issues for elu- 
cidating the superconducting properties. This problem 
is complicated by the structure of multiple Fermi sur- 
faces in LaOi-^F^FeAs systems, including hole Fermi 
surfaces otxp around the T-point at (k x ,k y ) = (0,0), 
and the electron Fermi surfaces fii^ around the M 12 
points at (7r,0) and (0,7r), respectively. Many theoretical 
proposals have suggested the fully-gapped extended s±- 
wave state which preserves the 4-fold rotational symme- 
try @, [3, [13, EH, E3, EDI- Experimentally, the superfluid 
density obtained from the penetration depth measure- 
ments is insensitive to temperature, consistent with this 
picture 14, [lj| Ell- Another competing pairing struc- 
ture in the square lattice is d x 2_ y 2 as proposed by sev- 
eral groups [3, 13, 17, HH. Kuroki et al. [3] find that 



the d x 2_ y 2 pairing wins over the s± state if the pairing 
contribution from a-Fermi surfaces is suppressed either 
by lowering the height of the As atom to the Fe plane 
or by heavily electron doping [I3| . In particular, several 
theoretical studies have sug gest ed that s± and d x 2_ y 2 
pairings are nearly degnerate[19l I20L l2l| . 

The calculation of the superconducting susceptibil- 
ity in the five-band model shows that the extended s- 
wave and d x 2_ y 2 pairing orders compete with each other 
[3|. The spin susceptibility has peaked values around 
(k x , k y ) = (tt, 0) and (0, 7r) which corresponds to the nest- 
ing wavectors connecting a and (3 Fermi surfaces, and 
also at nesting wavevectors around (n, tt/2) and (tt/2,tt) 
connecting two /3-surfaces. In the spin fluctuation mech- 
anism, the pairing order parameters favor opposite signs 
on two Fermi surfaces connected by nesting wavevectors. 
The first nesting favors the extended s-wave pairing, with 
the opposite signs for the electron and the hole fermi 
pockets, and the second one favors the d x 2_ y 2-wave pair- 
ing, with the opposite signs for the nearest-neighbor elec- 
tron pockets, as depicted in Fig. [TJ As one can see di- 
rectly from Fig. [TJ the pairing interactions based on these 
two different nesting vectors lead to a frustrating of the 
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FIG. 1: The two competing pairing order parameters 
LaFeAsOi-zFa, systems investigated in Ref. [8|, |9| which are 
compatible with the square lattice geometry: A) the extended 
s-wave (s±) and B) the d-wave (d x 2_ y 2). The nesting vectors 
are indicated by the bi-directional arrows. 



pairing order parameters - the pure extended s wave and 
the d x 2_ y 2 pairing states can not satisfy both nesting 
vectors simultaneously. In this situation, there arises a 
natural possibility of a mixed s + id pairing state, which 
can strike a compromise between the two nesting vectors. 

In this article, we investigate the possibility 
of time-reversal (TR) symmetry breaking states in 
LaOi-^F^FeAs systems based on the competition be- 
tween the extended s± and d x 2_ y 2 order parameters. 
Based on a Ginzburg-Landau (GL) free energy analy- 
sis, we show that after the occurrence of the s±-pairing 
at T c , the s + id pairing can develop at a lower temper- 
ature T' by breaking spatial rotation and TR symmetry 
under certain conditions. Spatial inhomogeneity can gen- 
erate supercurrent around non-magnetic impurities and 
the corners of square samples due to a symmetry allowed 
quadratic gradient coupling of these two order param- 
eters. The corner tunneling Josephson junction is an- 
alyzed. The resonance mode connecting two pairing or- 
der parameters can be measured through the B\ g Raman 
spectroscopy. 

Unconventional superconductivity with TR symme- 
try breaking effect has been actively investigated in re- 
cent years. Much experimental evidence has shown that 
Si'2Ru04 exhibits TR symmetry breaking effects consis- 
tent with the y x ± ip y state, including the /^SR and the 
Kerr effect 22, 13] • Recently, neutron scattering exper- 
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iments indicate the existence of the loop current in the 
pseudogap region 24, 25, [26| , and the Kerr effect has 
also been observed in the YBa 2 Cu30 6+2; system below 
the pseudogap temperature [27| . It has been pointed 
out that in multi-band superconductors TR symmetry 
breaking effect can exist with conventional pairing mech- 
anism due to the repulsive interband Cooper pair scatter- 
ings 28, 2ti 13, [3l|. LaOi-^F-rFeAs systems naturally 
have multi-band structure, thus it would be interesting 
to investigate possible consequences of the TR symmetry 
breaking pairing states. 

We construct a GL equation to describe the competi- 
tion between these two singlet pair channels of the ex- 
tended s-wave (s±) and <i-wave (d x 2_ y 2) order parame- 
ters as 



t =0.1 



F 



= a s (T)A* s A s + a d (T)A* d A d + Pi\A s \ 4 

+ ,3 2 |A d | 4 + /3 3 |A s | 2 |A d | 2 + f3 i ^A*A*A d A d + c.c.} 

= a s (T)A* s A s + a d {T)A d A d + &|A S . | 4 + (3 2 \A d \ 4 
+ ((3 3 + 2(3 4 )\A s ,\ 2 \A d \ 2 -(3 4 L 2 , (1) 



where a s « N \n(T/T s ) and a d w N ln(T/T d ); N is 
the density of states at the Fermi energy; (3s are not 
sensitive to temperatures; L = i(A*A d — A d A s ). Be- 
cause A s and A^ belong to different representations of 
the lattice symmetry group, they do not couple at the 
quadratic level but can couple at the quartic level. In 
particular, the /34-term is allowed by symmetry because 
orbital angular momentum is conserved modulo 4. The 
thermodynamic stability condition requires that j3\ t 2 > 
and 4/3i/3 2 - (/3 3 - 2|/3 4 |) 2 > 0. p 3 is also expected to 
be positive because of the competing nature between A s 
and A d . A similiar GL theory has been discussed in the 
high-T c cuprates(3^|. 

We assume that the pairing tendency in the s± is 
stronger but very close to that in the d x 2_ y 2 channel, 
i.e., T s > T d and 1 — T d /T s <C 1. By minimizing the GL 
free energy of Eq. [I] the condition for the appearance 
of Ad at temperature T is: —a d (T) > —a s (T)X, where 
^ = {k~ 2|/?4|)/(2/?i). This condition can be satisfied at 
A < I which gives rise to another critical temperature T" 
as: T = ^/T d T^(T d /T s )^^ < T d below which A d de- 
velops. The sign of p± determines whether TR symmetry 
is broken or not. If /3 4 > 0, it favors a phase difference 
of ±? between A s and A^, i.e., the s + id pairing. On 
the other hand, the real combination of s + d is realized 
at Pi < 0, which preserves TR symmetry. In both cases, 
the 4-fold rotational symmetry is broken at T < T' which 
corresponds to an Ising transition. In the former case a 
combined rotation of 90° and TR operation still leave the 
system invariant. 

While the uniform components of the A s and the A d 
components do not couple at the quadratic level, their 
gradient terms can. The general gradient terms of the 
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FIG. 2: Supercurrent induced by the impurity at f = (0, 0) 
for 04 > 0. We introduce the length scale £ = */ ls/\ot s | , 

energy scale Ao = ^\a s (T)\/2f3i, and the density of states 
No. The dimensionless parameters are defined as r'i = 

= ati/No, f3[ = ftAg/iVb, ^ = ~fi/£, 2 N Q , with the values 
used here given by r' = 0.1, a' = 10, a' s (T) = -1.0, a' d (T) = 
-0.75, p{ = 1.0, & = 1.0, p' 3 = 0.5, P' A = 0.4, 7 ^ = 1.0, j' d = 
0.5, and ^' ad = 0.35. The supercurrent is plotted in unit of 
jo = ea2 s(T)£,/ Pi and the length of each arrow is proportional 
to the magnitude of the supercurrent. and Cg> indicate the 
vorticities of the loop currents appearing in the four quadrants 
near the impurity. 



GL free energy are given by 

F g rad = 7s |nA s | 2 + 7 d|nA d | 2 + 7s(i (n;A:n ;c A 



- n; A*n y A d + c.c), 



(2) 



where II = V — 2ieA and A is the magnetic vector po- 
tential; the 7 S[ ;-term describes the coupling between the 
A Sj d orders allowed by the square lattice structure [33| . 
Minimizing the free energy we arrive the coupled equa- 
tions of 



a s A s + 2/3x|A s | 2 A s + p 3 \A d \ 2 A s + 2p 4 A 2 A* s 

- 7s (n 2 + n 2 )A s - lsd (ul - n 2 )A d = o, 

a d A d + 2/3 2 |A d | 2 A d + p 3 \A s \ 2 A d + 2/3 4 A 2 A* 

- 7d (n 2 + n 2 )A d - 7sd (n 2 - n 2 )A s = o. (3) 

The electric supercurrent can be represented as 

SF -. 
J = =]s+3d+3ds 
SA 

j s = 2ie ls [A* s UA s }+c.c, 
j d = 2ie ld [A* d UA d ] + c.c, 

Ids = 2ie lsd {(Ain x A d + A* d Il x A s )e x -(A* s Il y A d 
+ A* d U y A s )e y }+ c.c, (4) 

where j Std are the intra-component supercurrent, and j sd 
is the inter-component supercurrent. The total supercur- 
rent j satisfies the continuity condition: V • j = 0. 
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FIG. 3: Supercurrent distribution of a square sample in terms 
of the unit of jo with the same parameters as in Fig. [2] 
Current loops develop with the positive chirality at the right- 
down and left-up corners and the negative chirality at the 
other two corners. 




FIG. 4: Supercurrent components j s , jd and j s d in a sample 
of half plane geometry (x > 0) plotted with the unit of jo. 
The total current j — j s + jd + jsd vanishes. The parameters 
are the same as in Fig(2] 



A novel consequence of the 7 S( j-term is that spatial in- 
homogeneity induces supercurrent in a s + id supercon- 
ductor. We consider a non-magnetic impurity located 
at fo = (0, 0) modeled as a Gaussian distribution of as 
a s (r,T) = a s , d (T)+a e- r2 / r « with a > 0. At T < T', a 
s+id superconductor has the relative phase 9 s d = ±f be- 
tween A s ^ in the homogeneous system. However, in the 
presence of impurities, the spatial inhomogeneous distri- 
bution of order parameters causes the inhomogeneity of 
the relative phase 9 s d due to the j s d term, which induces 
supercurrent near the impurity. We solve Eqs. [3] and 
plot the supercurrent pattern in Fig. [2j Assuming the 
magnetic field generated by the impurity-induced super- 
current is small, we neglect the dependence of j on the 
magnetic vector potential A. The magnitude of the su- 
percurrent rapidly decays beyond the order of the healing 
length defined £ = . / -^f-,. The suppercurrent pattern is 



symmetric under the rotation of 180°, and a combined 
operation of TR and the rotation of ±90°. Along the x 
and y-axis passing the impurity, the supercurrents flow 
along these axis and exhibit the pattern of "two in" and 
"two out" , which is consistent with the reflection symme- 
try respect to the x and y-axis and the continuity condi- 
tion. In addition, the system has the combined TR and 
reflection symmetry respect to the diagonal axes of ±45°, 
thus supercurrents flow perpendicularly to these axes. 
This results in staggered current loops and thus small 
staggered fluxes in the four quadrants near the impurity. 
However, these current loops arise from both amplitude 
and phase variations and do not possess quantized fluxes, 
and thus are not of topological nature. At T c > T > T', 
a similar effect also occurs since a non-zero ci-wave or- 
der parameter can be induced despite the supercurrent 
in this temperature range is very small. On the other 
hand, if a s + d mixing is realized below T" at (3^ < 0, su- 
percurrents will not be induced by spatial inhomogeneity 
because s + d does not breaking TR symmetry. The spa- 
tial inhomogeneity only causes the amplitude variation 
of A Sl rf, whose relative phase difference remains fixed at 
sd = or 7T. 

We next consider a square sample and investigate the 
supercurrent at the boundary. The spatial variation of 
A s ,d is along directions normal to edges. The continuity 
condition suppresses the supercurrent except at the four 
corners as depicted in Fig. [3l Each of four corners de- 
velop a circulating supercurrent loop whose chirality are 
staggered as we move around the edges. This is also con- 
sistent with the combined symmetry operation of TR and 
the rotation of 90°. Because at each corner the supercur- 
rent has the same chirality, thus is easier to be detected 
by using SQUID than the single impurity case. 

Furthermore, we consider the spatial distribution of 
A Sy d at 04 > along the edge of a half-plane sample 
with x > and the boundary of the y-axis. Given the 
boundary conditions of A s (x = 0, y) — Ad(x = 0, y) = 
and the spatial homogeneity along the y-direction, Eq. 
03 reduces to coupled ID equations. Again the conti- 
nuity condition forbids the appearances of a net super- 
current for a s + id superconductor, although the three 
components of the supercurrent, j s , jd, and j s d, are non- 
zero individually as depicted in Fig. [4] One important 
feature is that j s and jd have opposite signs. This in- 
ternal counterflow supercurrents is a result of the ^ s d 
term. By plugging A St d(x) = p Sy d{x)e tes - d ^ into the 
GL free energy, such a term can be approximated as: 
2j sd sin(9 s - d )[(d x 9 s )(d x p d )ps - (d x d )(d x p s )p d }, and 
therefore the free energy can be lowered by choosing 
sgn(d x 9 s ) = -sgn(d x 6 d ). 

The Josephson junctions touching edges with different 
orientations of a single crystal sample act as convinc- 
ing phase sensitive measurements for the ci-wave pairing 
symmetry in high T c systems I34l . |35 | and p-wave pair- 
ing symmetry in Sr2Ru04 36|, 37j . We consider the same 
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cornier junction in which the two adjacent faces of a single 
crystal FeAs superconductor are connected via Josephson 
weak coupling with a conventional s-wave superconduct- 
ing thin film. We define that J Sy d are the Josephson cou- 
pling between the s± (d) order parameter and the conven- 
tional s-wave order parameter, respectively. Assuming 
that the two junctions at two adjacent faces are equal, 
we express the critical current with applied magnetic 
flux § ext asjUj]: I c (<f>ext) = 2I Q cos[n ($ ext /$ ) + 5/2], 
where 6 = tan -1 [2ej sin 9 s d/(l — e 2 )] , ej = Jd/J s , and 
$o is the flux quantum. 5 could only be or ir for the 
s + d mixing depending on which component is dominate, 
while it can be any value between and it for the case of 
s + id mixing. 

Now we consider the case of (A s ) = Ao and (Ad) = 
0. In this case we predict that a new collective mode 
can be observed in the £?i ff -mode of Raman spec- 
troscopy, which behaves as the resonance mode con- 
necting A Sy d, and is approximately the nematic op- 
erator Nd associated with the B\ g mode in the Ra- 
man spectroscopy. One can check the following com- 
mutator: [N d ,A d ] = -2A\ + A|', where N d = 
l/VJ2 k ( coak x - cosk y )ctcr, A f d = l/VJ2k(cosk x - 



5 ki, ) c^— 

v > k] 



-H 



At 



i/v J2k cos k * cos fc «4 T -fci 



COS" 



k y^ c h c -ki- 



and 

Since (A s > ^ 



0, Ad and N d are conjugate variables which lead to collec- 
tive modes similiar to the n resonance mode in cuprates 

where K Nd ~ AqA 2 , and 



0J r 



V K N d 



(J 

K Ad ~ a* d (T) = a d + (/3 3 - 2|/3 4 |)Ag are the stiffnesses 
for N d and A^ respectively, and C ~ Nq is the Berry 
curvature between N d and A d - Following the same ar- 
guement in Ref. [H, |4lJ , the excited state with energy 
uj res can be defined as: \Ad) — A^|0), where |0) is the 
BCS ground state. As a result, at T — the Raman 
response function of the B\ g mode becomes: 



u) + uj n — ie _ 
|(0|[iV d ,At]|0)| 2 |(0|[iV d ,At]|0)| 2 ' 



IC 



IC 



(5) 



which has a sharp peak at uj = uj res . The temperature 
dependence of uj res may be complex. Nevertheless, the 
revelence of A^ can still be infered from it. If the A^ 
is competing to A s but the mixed state is not favored 
(A > 1), uj res should have weak temperature dependence 
and remain finite as T — > 0. On the other hand, if the 
mixed state could occur at T = T' (A < 1), u res should 
appear after T < T c and decrease dramatically to zero as 
T is approaching T' . 

In conclusion, we have investigated the competition 
between the extended s-wave and d-wave Cooper pairing 



orders in the FeAs-based superconductors. The multiple 
nesting wavevectors naturally leads to the possibility of a 
s+id pairing state which breaks the TR symmetry. Based 
on a general Ginzburg-Landau theory we have shown that 
in such a state supercurrent can be induced by spatial 
inhomogeneity, and several possible experiments to de- 
tect this state are discussed. We also proposed that a 
new collective mode should be observed in the B\ g Ra- 
man spectroscopy as the resonant mode between the two 
competing order parameters. 
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